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ABSTRACT 



Galaxy-galaxy interactions are expected to be responsible for triggering massive star formation and possibly accretion onto a super- 
massive black hole, by providing large amounts of dense molecular gas down to the central kiloparsec region. Several scenarios to 
drive the gas further down to the central ~ 100 pc, have been proposed, including the formation of a nuclear disk around the black hole, 
where massive stars would produce supernovae. Here, we probe the radial distribution of supernovae and supernova remnants in the 
nuclear regions of the starburst galaxies M82, Arp 299- A, and Arp 220, by using high-angular resolution (< 0."1) radio observations 
published in the literature (for M82 and Arp 220), or obtained by ourselves from the European VLBI Network (Arp 299-A). Our main 
goal was to characterize the nuclear starbursts in those galaxies and thus test scenarios that propose that nuclear disks of sizes ~100 
pc form in the central regions of starburst galaxies. We obtained the radial distribution of supernovae (SNe) in the nuclear starbursts 
of M82, Arp 299-A, and Arp 220, and derived scale-length values for the putative nuclear disks powering the bursts in those central 
regions. The scale lengths for the (exponential) disks range from ~20~30 pc for Arp 299-A and Arp 220, up to ~ 140 pc for M82. The 
radial distribution of SNe for the nuclear disks in Arp 299-A and Arp 220 is also consistent with a power-law surface density profile 
of exponent y = 1, as expected from detailed hydrodynamical simulations of nuclear disks. Our results support scenarios where a 
nuclear disk of size ~ 100 pc is formed in (U)LIRGs, and sustained by gas pressure, in which case the accretion onto the black hole 
could be lowered by supernova feedback. 
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1. Introduction 

The activity in the central regions of luminous and ultra- 
luminous infrared galaxies (LlRGs, IO'^Lq ^ ^ir < lO'^L©, and 
ULIRGs, LiR > lO'^Lo, respectively) is powered by either accre- 
tion onto a supermassive black hole (SMBH), and/or by a burst 
of activity due to young, massive stars. It has been found that 
LlRGs and ULIRGs are associated with interacting and merging 
galaxies, respectively, and the interaction between the galaxies 
is assumed to provide large amounts of dense molecular gas that 
eventually reach the central regions of the galaxies, where they 
are responsible for triggering mas sive star formation and possi- 
bly accretion onto an SMBH (e.g. JPi Matteo et al. 2007) . While 
interactions and bars seem to aid the transport of the molecu- 
lar gas into the kiloparsec region of galaxies (e.g., 'Simk in et al.l 
[l980), it has been recognized that a major difficulty is in driving 
the gas from the kiloparsec region down to the central parsec, 
or parsecs, of (U)LlRGs, where a nuclear starburst and/or an ac- 
tive galactic nucleus (AGN) are expected to exist. Several mech- 
anisms have been proposed to overcome this problem, includ- 
ing tidal torques driven by either ste llar bars (Shlosman et al. 



199C) or major and minor mergers (Mihos & Hernquistlll994[ 



1996 *). gas drag and dynamical friction in a dense n uclear star 



clust er (Norman & Scoville 1988), or radiation drag (lUmemural 
l2001tlKawakatu & Umemurall2002l) . However, not all the accu- 
mulated gas is expected to be accreted directly onto the puta- 
tive SMBH, since the removal of the a ngular momentum of the 
gas is very inefficient (IUmemural200 II) . The angular momentum 



that is not removed will then permit the accreted gas to form 
a reservoir, na mely a nuclear disk in the central ~100 pc around 
the AGN (e.g. JWada & Normanl2002llKawakatu & Wadal2008l 
and references therein). Furthermore, if the gas in the nuclear 
starburst is dense enough, vigorous star formation is expected 
to occur, which can be detected via, e.g., high-angular resolu- 
tion observatio ns of recently exploded supernovae (e.g., in the 
ORG Arp 299 dPerez-Torres et al.l2009l) . or the ULIRG Aip 220 
(iLonsdale et al.ll2006l) '). The nuclear starburst is also expected to 
affect the growth of the SMBH, since the radiation and/or super- 
nova (SN) feedback caused by starbursts can eventually trigger 
the mass accretion ont o a SMBH (e.g.. lUmemuraet al]|1997t 
IWada & Normanll2002h . 

iKawakatu & Wadal ( l2008l) proposed an scenario for 
(U)LIRGs, where a nuclear disk is supported by the turbulent 
pressure of SNe. In their model, the turbulence due to supernova 
activity transports the angular momentum toward the central 
pc -scale region. One o f the implications of the model developed 
bv lKawakatu & Wadal ([2008) is that the BH growth rate depends 
on the spatial distribution of both young stars and SNe in the 
nuclear disk. However, the study of this spatial distribution is 
a challenging observational task for two main reasons. First, 
the huge column densities present in nuclear starbursts cause 
enormous extinction towards the central regions of (U)LlRGs 
at optical and infrared wavelengths, such that essentially no 
supernova would be detected at those wavelengths in the central 
few hundred parsecs of a (U)LIRG. Second, the enormous 
angular resolution needed to discern individual objects at the 
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typical distances to (U)LIRGs ( > 50 Mpc) requires milliarc- 
second angular resolution to resolve individual core-collapse 
supernovae (CCSNe) and/or supernova remnants (SNRs). 

In this paper, we use high-angular radio observations from 
the literature, as well as data obtained by us, of three starburst 
galaxies in the local universe (M82, Arp 299-A, and Arp 220) to 
analyze the radial distribution of CCSNe/SNRs in their central 
few hundred parsecs, which may have strong implications for 
the (co)-evolution of AGN and (nuclear) starbursts in galaxies, 
in general, and in (U)LIRGs in particular. 



2. The radial distribution of SNe in nuclear 
starbursts 

The radial distribution of SNaj on galactic scales has been thor- 
oughly studied (see, e.g., Bartunov et al. 1992; van den Berghl 
[1997; Anderson & James 2009; Hakobyan et al. 2009). Optical 
observations have been very useful for this purpose, as dust 
extinction plays a minor role in the external regions of most 
galaxies. However, the dust-enshrouded environments of the 
nuclear regions in (U)LIRGs, as well as the need for mil- 
liarcsecond angular resolution to pinpoint individual SNe, have 
prevented this kind of studies in the central regions of local 
(U)LIRGs. Fortunately, both issues -obscuration and (angular) 
resolution- can be overcome by means of Very Long Baseline 
Interferometry (VLBI) radio observations of the central regions 
of (U)LIRGs, since radio is dust-extinction free, and VLBI pro- 
vides milliarsecond angular resolution. 

2.1. Methods 

We used VLBI observations to probe the radial distribution 
of SNe in the nuclear regions of Arp 299-A (3 1 objects) and 
Arp 220 (48 objects), as well as MERLIN, VLA, and VLBI ob- 
servations towards M82 (38 objects). For each galaxy, we fit- 
ted the surface density profile of the SNe to two different disk 



profiles: (i) an exponential disk, E 



S^™ exp (-r//isN), and 



(ii) a disk with a power-law density profile with radius, I,^^ - 
^o'^('"/'"out) ''■ The first case (exponential disk) assumes that the 
SN distributi on follow s the radial surface-brightness density in 
spiral disks ( Freemanl ll970). while the second one (power-law 
disk) was used to probe the profiles predicted by numerical sim- 
ulations (Wada & Norman 2002). The notation follows the one 
used in|Hakobvan et al. (2009), where Eg'^ is the surface density 
of SN at the center, r is the normalized radius, and rout is the ra- 
dius of the farthest SN analyzed, considered here to be the outer 
boundary of the putative nuclear disk. The main parameters that 
are obtained from those fits are the scale length, /zsn, and the in- 
dex of the power-law profile, y, of the putative nuclear disks in 
the central regions of these galaxies. 

Our approach consisted in model fitting the surface density 
profile of the radial distribution of SNe, by determining the num- 
ber of sources in concentric rings around the center of the galaxy. 
For this purpose, we carried out a simil ar data analysis to the 
one described in iHakobvan et al.l (l2009l) (hereafter, H09), and 
we refer the reader to that paper for further details. In partic- 
ular, we calculated the de-projected distance of each compact 
source to the center of the galaxy (the AGN, whenever its iden- 
tification was possible) from the position angle of the galaxy 
major axis, and the inclination of the galactic disk. We obtained 
those parameter values from the HyperLeda database. We also 



followed H09, and measured radial distances in units of R25, 
i.e., the isophotal radius at which the surface brightness was 
25 mag/arcsec" in the Johnson B filter We then derived the sur- 
face density profile by determining the number of sources in con- 
centric rings around the nucleus, and fitting the data to both an 
exponential and a power-law disk, as described above. 

Since the supernova samples for all our three galaxies are of 
relatively small size, resulting in some small size bins, the actual 
uncertainties in the fits could be larger than formally obtained 
from a single fit. To assess the robustness of our fitting proce- 
dure, we performed a series of Monte Carlo (MC) simulations. 
For each galaxy, we generated 10000 mock supernova samples, 
whose values were uniformly distributed within the uncertainty 
of each data point, and fitted them. In this way, we obtained a 
distribution of 10000 different fitting values for the scale-length 
parameter, /isN, and the power-law index, -y, for each galaxy. This 
approach allowed us to obtain more reliable values of those pa- 
rameters, since no assumption was made about the Gaussianity 
of the distributions of both the scale length and power-law in- 
dex, or of their errors. We then used the median value of the 
distribution of our MC values to characterize /isN and 7, and 
set the uncertainty in these parameters at the 90% confidence 
level. To test our method, we reanalyzed the sample of H09 for 
239 CCSNe within 216 host galaxies, and obtained a value of 
/zsN - 0-29^oyj, which is in excellent agreement with their esti- 
matqjl and confirms the robustness of our analysis. 



2.2. Results 

We outline the results obtained for each individual galaxy. We 
refer the reader to the Appendices for details on the SN data 
used in our fitting analysis. 



2.2.1. Arp 299-A 

Arp 299-A (D - 44.8 Mpc) is one of the two interacting galax- 
ies of the LIRG Arp 299. The sample of SNe used for our 
study comes from both observations from t he European VLBI 
Network (E VN) at 5.0 GHz, p resented in iPerez-Torres et al.l 
(2001) and in lBondietal] ( 1201 21) (see also Appendix|Al), as well 
as from published data obtained with the Very Long Baseline 
Array (VLBA) at 2.3 and 8.4 GHz bv Ulvestad (2009) . The ex- 
istence and precise loc ation of the low luminosit y AGN in Arp 
299-A was reported bv lPerez-Torres et al.l (l20I0l) . hence its po- 
sition was used to calculate accurate de-projected distances for 
all SNe in the nuclear region of Arp 299-A. The sample has a 
total of 31 sources. However, although there are sources at dis- 
tances as far as r = 1 1 .4 x 10"^, only one source is further than 
f = 6.6 X 10"^, so the latter value of ? was used as the limiting 
radius for our study, giving a final number of 30 sources. The 
scale length obtained using the non-linear fit is /isn = 1 .6 x 10""', 
with a standard deviation of 0.5 x 10"^. After performing the 
MC analysis, we obtained a final value /jsn of (l-9;^Qg) x lO"-' 
and S^^ = (1.7+2-8) xlO^ 

2.2.2. Arp 220 

Arp 220 is a merging system at 77 Mpc, and the closest and 
most well-studied ULIRG. High-angular resolution radio ob- 



In this paper, the term SNe denotes both CCSNe and SNRs. 



- H09 report /isn = 0.29 ± 0.01. We note that the explicit list of 
CCSNe is not given in H09, and therefore we could not reproduce ex- 
actly their sample. 
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servations from its tw o nuclear regions were first obtained by 
iLonsdale et alJ (|2006) using global VLBI, who found a total of 
49 radio supernovae in both the east (20) and west (29) nuclei. 
We used the position of the peak flux density from a delay-rate 
VLBI map as the estimated p osition for the center of each galac- 
tic nucleus (iParra et al.ll2007l) . One of the compact sources in the 
Arp 220 east nucleus is at a much larger distance than the rest 
of them, which we did not use in our fit, thus limiting the nor- 
malized radii to r = 0.011. Therefore, a total of 48 SNe were 
considered in the fits, which we carried both separately, for each 
nucleus, and for the combined data set (Fig. IA.3b . Similarly, 
we performed MC simulations for each of the two nuclei, as 
well as for the combined sample In this way, we obtained a 
scale-length for the east and west nuclei of (S.l^y^j x 10"^ and 
(3-4*Qg) X 10"^, respectively, and a combined scale-length for 

Arp 220 is (3.3!{}^) x 10^1 To quantify the effect of the uncer- 
tain position of the putative AGNs, we shifted them 100 times 
randomly around our nominal position, given by the peak flux 
density, by a maximum value of o-^kj^dec) - FWHM/(2 S/N), 
where FWHM is the full-width at half maximum, and S/N is 
the signal-to-noise ratio, resulting in 4.2 mas for the east nu- 
cleus and 1 .4 mas for the west one, and then fitted to the data. 
The maximum differences obtained were considered as a sys- 
tematic error, which we added to our nominal values, such that 
the combined scale length became Jisn - w-^lol) ^ lO""', with 



Arp 299-A and Arp 220 E+W combined 



ySN 
^0 



(7.6!|?) X 10\ 



2.2.3. Combination of Arp 299-A and Arp 220 

The normalized radii of Arp 299-A and Arp 220 (both east and 
west nuclei) are comparable. We can thus combine every source 
prior to the analysis to obtain an average scale-length parameter 
with smaller errors to help characterize both nuclear starbursts. 
The scale length obtained for the combined samples is /isN - 
(2.0^Q^)xlO"^. The original non-linear fit is shown in Fig. [Tjand 
our MC analysis results are shown in Fig. |2] After taking into 
account systematic errors due to the uncertainty in the precise 
location of the AGNs of the Arp 220 east and west nuclei, the 
final value for the scale length is /zsn - {^■^^qI) ^ 10"^^ and 

ySN 
^0 



(3.0!°:?) X 10^ 




0.004 0.006 0.008 0.010 



Fig. 1. Non-linear fit of the surface density profile of the com- 
bined sample of SNe from both Arp 299-A and Arp 220. 



Arp 299-A and Arp 220 E+W combined 




0.0014 0.0015 0.0018 0.0020 0.0022 0.0024 

Fitting value 



0.0025 0.0028 



Fig. 2. MC histogram of /isN for the combined sample of SNe 
from Arp 299-A and the two nuclei of Arp 220. The original 
non-linear fit is shown in red, and the median of the distribution 
created with the MC simulation is shown in green. The dashed 
lines correspond to the 90% confidence level interval. 



2.2.4. M82 

At a distance of 3.2 Mpc, M82 is the prototypical starburst 
galaxy, which has been the target of many observations at 
high-angular resolution. For our analysis, w e combined source s 

(Il997h 



WiUs et al. 



catalogued as SN e (or S NR) detect ed b\ 

[Allen & Kron berg(IT998l). [M cDonald et al.1 (l2002h . lFenech etd] 
O2008i) . and iBrunthaler et al.i (i2009). These comprise observa- 
tions at several frequencies with MERLIN, VLA, and VLBA, 
including a total of 39 sources. Since no AGN has yet been 
discovered i n M82, we instead use d the radio kinematic cen- 
ter of M82 (IWeliachew et all [19841) as the position of its nu- 
cleus. The scale length obtained with the MC simulation is 
/isN = (^-S^os) ^ 10"^. Although the radio kinematic center is 
likely to be very close, if not exactly coincident, with the puta- 
tive nucleus of M82, we nevertheless quantified the effect that an 
error in the position chosen as the nucleus of M82 would have in 
the fits. In particular, we shifted our fiducial position for the nu- 
cleus by up to 20 parsec in several directions, and then fitted the 



data. The maximum differences were considered as a system- 
atic error, which we added to our nominal value, thus yielding 
/2SN = (2.8!°:^) X 10"'' and S^^ = (8.3!|^) x 10^ 

3. Discussion and summary 

Our main results are summarized in Table[T]and Fig. [3] where we 
show the scale lengths, index of the power-law SN profile, and 
the radial distribution of SNe for the three galaxies discussed in 
this paper, as well as for the sample of spiral galaxies from H09. 
We note that the scale lengths for the exponential disk fits to 
the radial distribution of SNe in the nuclear regions of the star- 
burst galaxies, are at least an order of magnitude smaller that 
those obtained by H09 for the whole disks of spiral galaxies. In 
particular, /zsn is two orders of magnitude smaller in the case of 
Arp 299-A and Arp 220. Correspondingly, the physical sizes in- 
ferred for the scale lengths of the nuclear disks (col. 3 in Table[T]) 
are much smaller than for galactic disks: ~20-30 pc for Arp 299- 
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Table 1. Scale length parameters for the galaxies discussed in 
this paper 



Nucleus 



/isn/IO-' 



^SN (pc) 



Arp 299-A 


1 Q+1.9 




9Q Q+29.6 


1 r°2 


Arp 220 East 


3 1+2° 

■^•^-0.9 




99 9+14.4 
^^•^-6.2 


1.0!- 


Arp 220 West 


3 4+'« 




74 4+11.2 
^^•^-10.8 


OKI 


Arp 220 E+W 


^■ni 




23.4!^:^ 


o.8!°i 


A299 + A220 


2-0!gi 




20.9!?* 


o.9!°:l 


M82 


(2.8!-) X 


10' 


1 44.4^21 -5 


- 


H09 sample 


(2.9-?) X 


10^ 


- 


- 



All quoted uncertainties correspond to 90% confidence limit inter- 
vals, to which a systematic error was added in the cases of Arp 
220 east, Arp 220 west, and M 82, owing to the uncertainty in the 
precise location of their centers. 



A or Arp 220, which is similar to the size derive d for the nu- 
clear s tarburst of the LIRG NGC 7469 (30-60 pc: lDavies eFaP 
l2004bl) . using CO interferometric observations, and ~160 pc for 
M82, which is also similar to t he scale-length of th e nuclear disk 
in the ULIRG/QSO Mrk 231 (iDavies et al.ll2004ah . 

We show in Fig. |3]the radial distribution of SNe in our three 
starburst galaxies and, for comparison, also for spiral galaxies 
(from the H09 sample), along with the fits to these data. We 
note the clear existence of three different regimes: one charac- 
terized by a very steep profile of the surface number density 
of SNe, which is typical of strong, nuclear starbursts (Arp 299- 
A-nArp 220); a second, less steep profile, as indicated for the 
(circum)nuclear starburst in M82; and a third, flatter profile, 
which is typical of very large disks, such as those in spiral galax- 
ies (H09 sample). These results suggest that the surface density 
of SNe, hence of available gas to convert into stars, reaches a 
maximum in the vicinities of the SMBH in LIRGs and ULIRGs. 
We also show in col. 4 of Table [T] the fits to a power-law disk 
profile, consistent with the fiducial value of y = 1 used by 
fWada & Norman (2002) and iKawakatu & Wadal (l2008h for the 
nuclear disk, and appears to yield further support to their mod- 
eling. 

We have modeled for the first time the radial distribution of 
SNe in the nuclear starbursts of M82, Aip 299-A, and Arp 220, 
and derived scale-length values for exponential disks, which are 
in the range between ~20 pc and ~140 pc. We have also modeled 
these SNe distributions as power-law disk profiles, and found 
that they agree with state-of-the art numerical simulations of nu- 
clear disks around SMBHs. We interpret our results as evidence 
of nuclear disks around the centers, i.e., AGNs, of starburst- 
dominated galaxies and also support for evolutionary scenar- 
ios where a nuclear disk of size < 100 pc is formed in LIRGs 
and ULIRGs ( Kawakatu & Wada 200 8). In particular, the scale- 
length obtain ed for the LIRG Arp 299- A and the LLAGN nature 
of its SMBH dPerez- Torres et al.ll2010h . suggests that its nuclear 
disk is likely supported by gas pressure, such that the accretion 
onto the SMBH is smaller than for a turbulent-supported disk. It 
is very likely that this is also the case for the ULIRG Arp 220. 
Future deep VLBI observations of a significantly large sample 
of (U)LIRGs, which would result in the discovery of new SN 
factories, will be of much use for deriving statistically signifi- 
cant results on the size of nuclear disks, and thus setting useful 
constraints on (co)-evolutionary scenarios. 




■ Arp 299 + Arp 220 

• M82 

♦ H09 sample 



0.05 0.10 0.15 0.20 0.25 

Fig. 3. Radial distribution of SNe for (nuclear) starbursts and spi- 
ral galaxies. The data correspond to the combined Arp 299-A -i- 
Arp 220 sample (blue squares), M82 (green circles), and the H09 
sample of spiral galaxies (red diamonds), and the fitted lines are 
for Arp 299-A + Arp 220 (soHd fine), M82 (dotted-dashed line), 
and H09 (dashed line). For the sake of clarity, only the innermost 
part of the H09 sample is shown in the plot. 
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Appendix A: Arp 299-A and Arp 220 



Table A.l. Arp 299-A" source list 



Arp 299-A - Power law fit 



Name 


A a" 


AS 


Deprojected dist. 




(J2000.0) 


(J2000.0) 


(mas) 


A3 


33.5941 


46.560 


252.3 


A15 


33.5991 


46.637 


190.9 


33.612+46.70'^ 


33.6124 


46.696 


71.3 


33.615+46.67'^ 


33.6154 


46.667 


47.9 


A26 


33.6163 


46.652 


55.6 


33.617+46.72'^ 


33.6173 


46.724 


44.2 


A18 


33.6195 


46.463 


276.6 


A19 


33.6198 


46.660 


45.1 


A6 


33.6207 


46.597 


122.6 


AO 


33.6212 


46.707 


13.2 


A2 


33.6218 


46.655 


60.5 


A27 


33.6221 


46.692 


25.06 


33.624+46.77^ 


33.6241 


46.771 


81.9 


33.627+46.44^ 


33.6267 


46.440 


331.2 


A28 


33.6272 


46.719 


66.8 


33.628+46.62^ 


33.6280 


46.623 


135.8 


33.629+46.65"^ 


33.6290 


46.647 


121.6 


A7 


33.6301 


46.786 


117.7 


AS 


33.6305 


46.401 


394.6 


A9 


33.6307 


46.620 


159.4 


A29 


33.6344 


46.915 


249.6 


A22 


33.6355 


46.677 


160.6 


A30 


33.6362 


46.863 


209.1 


AlO 


33.6392 


46.551 


291.6 


All 


33.6403 


46.582 


272.9 


A31 


33.6441 


46.629 


271.1 


A12 


33.6495 


46.590 


347.3 


A4 


33.6500 


46.537 


394.0 


A25 


33.6523 


46.701 


311.8 


A13 


33.6531 


46.733 


309.6 


A14 


33.6825 


46.570 


666.2 




Fig. A.2. Surface density of SNe vs radii normalized to the fur- 
thest source, rout, in Arp 299-A and non-linear fit to a power law 
with y = 1.1. 
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'' Arp 299-A has a major axis position angle of PA = 163.6° and an 

inclination of its galactic disk of ; = 52.4° . 
'' Coordinates are given with respect t o a(J2000.0) = 11 :28:00.0000 

and 5(J2000.0) = 58:33:00.000, (see lBondi etaDIM^ . 
■^ Sources from iUlvestad (i2009h . 



Arp 299-A 



Fig. A.3. Non-linear fits to the surface density profile of the ra- 
dial distribution of SNe in Arp 220. Top left: Fit to the 19 inner 
sources in the eastern nucleus (/isN = 2.7 x 10"-'). Top right: 
Fit to the 29 sources in the western nucleus (/isN = 3.5 x lO"-'). 
Bottom: Combined fit for the SNe in the east and west nuclei (48 
sources), yielding Ji^n = 3.3 x 10"^ and 2^^ = 6.0 x 10^ We 
note that some of the points are in rings where no sources were 
detected. While their corresponding surface density is zero, their 
upper limits were used in the fit. 
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Fig. A.l. Surface density profile of the 30 inner SNe in Arp 299- 
A. The solid green line is a non-linear fit to an exponential disk 
with /isN = 1 .6 X 10-3 and ^^ = 1 .6 x 10^ 



Ruben Herrero-IUana et al.: Evidence of nuclear disks in starburst galaxies from their radial distribution of supemovae, Online Material p 2 



Arp 220 E+W - Parsec scale 
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Fig.A.4. Surface density profile of the combined sources for 
the east and west nuclei of Arp 220, in parsecs. The non-linear 
fit yields a scale length of /zsn - 23.4 pc, with I,^^^ = 1.2 x 
10-2 pc-2. 

As p 299-A and Arp 220 E+W combined - Parsec sc ale 
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Fig.A.5. Surface density profile of the combined sources for 
Arp 299-A and Arp 220, in parsecs. The non-linear fit yields 
a scale-length of /zsN = 19.4 pc and 2^'^ = 3.0 x 10"^ pc^^ 



Arp 220 East - Power law fit Arp 220 West - Power law fit 
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Arp 220 E+W - Power law fit 
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'' Aip 220 has an inclination angle of 57.0°, with a major-axis posi- 
tion angle of 96.5°. 

'' Coordinates are given with respect to a(J2000.0) = 15:34:00.0000 
and 5(J2000.0) = 23:30:00.000. 
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Fig. A.6. Non-linear power-law fits of the surface density of SNe 
vs radius, normalized to the furthest source, rout, in Arp 220. Top 
left: East nucleus, with j - 1.0. Top right: West nucleus, with 
y - 0.7. Bottom: Combined analysis for the East and West nuclei 
(48 sources), yielding y - 0.7. Note that some of the points are in 
rings where no sources were detected. While their corresponding 
surface density is zero, their upper limits were used in the fit. 



Ruben Herrero-IUana et al.: Evidence of nuclear disks in starburst galaxies from their radial distribution of supemovae, Online Material p 3 



Arp 299-A and Arp 220 E+W combined - Power law fit 




Fig. A.7. Surface density of SNe vs radii normalized to the fur- 
thest source, rom, for the combined sample of SNe from Arp 299- 
A and the two nuclei of Arp 220. A non-linear fit to a power law 
with radius yields y = 0.9. 
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Appendix B: M82 



Table B.l. M82^ source list 



Name" 


A a" 


A^ 


Deprojected dist. 




(J2000.0) 


(J2000.0) 


(mas) 


39.10+57.3 


47.87 


43.74 


37.03 


39.40+56.2 


48.16 


42.53 


29.20 


39.64+53.3 


48.41 


39.68 


20.00 


39.77+56.9 


48.54 


43.36 


27.84 


40.32+55.2 


49.07 


41.53 


16.98 


40.61+56.3 


49.32 


42.19 


15.93 


40.68+55.1 


49.42 


41.42 


14.13 


41.30+59.6 


50.05 


45.92 


20.40 


41.95+57.5 


50.70 


44.00 


6.94 


42.61+60.7 


51.35 


46.89 


9.35 


42.66+51.6 


51.38 


47.76 


13.21 


42.67+55.6 


51.37 


41.78 


17.44 


42.67+56.3 


51.40 


42.65 


13.42 


42.80+61.2 


51.54 


47.53 


10.25 


43.18+58.2 


51.91 


44.58 


8.70 


43.31+59.2 


52.03 


45.42 


5.75 


43.40+62.6 


52.11 


48.74 


10.22 


43.55+60.0 


52.27 


46.31 


4.22 


43.72+62.6 


52.43 


48.79 


7.63 


43.82+62.8 


52.54 


48.98 


7.71 


44.01+59.6 


52.72 


45.77 


12.03 


44.08+63.1 


52.75 


49.46 


8.50 


44.28+59.3 


52.99 


45.53 


16.33 


44.34+57.8 


53.05 


43.94 


24.65 


44.40+61.8 


53.13 


47.87 


8.35 


44.51+58.2 


53.22 


44.36 


24.59 


44.89+61.2 


53.61 


47.35 


15.69 


45.17+61.2 


53.88 


47.41 


18.61 


45.24+65.2 


53.96 


51.36 


12.15 


45.42+67.4 


54.13 


53.53 


18.13 


45.52+64.7 


54.22 


50.94 


12.94 


45.75+65.3 


54.44 


51.43 


14.17 


45.79+64.0 


54.50 


50.22 


16.24 


45.89+63.8 


54.60 


49.98 


17.85 


46.52+63.9 


55.21 


50.05 


24.19 


46.56+73.8 


55.26 


59.90 


38.39 


46.70+67.1 


55.41 


53.16 


19.71 


47.37+68.0 


56.10 


54.00 


24.01 


SN2008iz 


51.55 


45.79 


2.41 



^ M82 has a major position angle of 67.5° and an inclination of 79.4° . 
'' Names from the literature are derived from B1950 coordinates. 
■^ Coordinates are given with respect to q-(J2000.0) = 09:55:00.00 and 
^(J2000.0) = 69:40:00.00. 



Appendix C: Nature of the radio sources 

The sources used in this paper are classified generically as 
supemovae (SNe), by which we mean core collapse super- 
novae (CCSNe) and supernova remnants (SNR). In the case 
of M82, the closest of the three galaxies d iscus sed here, we 
used t he source list in McDonald et al.l (l2002h and lFenech et alJ 
(l2008h . where the authors already discriminate SNe from HII 
regions based on their spectral shape, brightness temperature, 
and observed shell structure. In the cases of the more dis- 
tant galaxies, Aip 299 and Arp 220, we used available spec- 
tral information from the literature, and especially, the bright- 
ness temperatures inferred from VLBI observations, which for 
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Fig. B.l. Surface density profile for the radial distiibution of SNe 
in M82 (39 objects). The solid green line is a non-linear fit to an 
exponential disk with /zsn = 2.5 x 10"^ and 2^'^ = 7.5 x 10^ 



M82 - Power law fit 




Fig. B.2. Surface density of SNe vs radii normalized to the fur- 
thest source, rout and non-linear fit to a power law (nuclear) disk 
for M 82. The poor goodness of the fit does not seem to be con- 
sistent with a nuclear disk whose SN surface density follows a 
power law in the central few hundred pc. 



all of the sources c orrespond to non-thermal radio emitters (e.g., 
Perez- Torres et al.l (2009 ) for Arp 299 A, Parra et al. (2007) for 
Arp 220). 

Single or multiple ultra compact HII (UC HII) regions can 
be ruled out because their ionizing radiation is normally domi- 
nated by a single O star, whose maximum free-free thermal ra- 
dio luminosity is several orders of magnitude below the radio 
luminosity from any of the objects detected in either Arp 299- 
A or Arp 220. The expected thermal radio emission from super 
star clusters (SSCs), which can host up to a few 10^ stars, sim- 
ilarly cannot be responsible for the observed emission. At the 
distances of Arp 299 and Arp 220, the VLBI beam sizes cor- 
respond to ~ 1.5 pc, implying radial sizes of at most about 0.8 
pc. Young massive stellar clusters of this size have a mass of at 
most 5 X 10^ Mq (e.g. Portegies Zwart et al. 2010). Assuming 
a Kroupa initial mass function with a = -1.0 between 0.1 and 
0.5 Mq, and a = -2.35 between 0.5 and 100 Mq, the number 
of 06 and earlier type stars is 142 (out of a total of 10^ stars), 
hence the maximum attainable thermal radio luminosity would 
be ~ 2 X lO^'' erg/s/Hz, which is at least an order of magnitude 
lower than observed in Aip 299-A or Arp 220. 



